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The set up of matched filters for the detection of gravitational waves from in-spiraling compact 
binaries is usually carried out using the restricted post-Newtonian approximation: the filter 
qq ■ phase is modelled including post-Newtonian corrections, whereas the amplitude is retained at 

(^-^ ' the Newtonian order. Here we investigate the effects of the introduction of post-Newtonian 

corrections also to the amplitude and we discuss some of the implications for signal detection 
If) • and parameter estimation. 

o 
o 
o 



cr 



X 



1 Introduction 



It has long been recognized that binary systems of compact objects are important sources of 
gravitational waves (GW) both for the ground based interferometric detectors currently under 
construction &@ and future space-based interferometers like LISAeI. 

In this paper, we consider the in-spiral phase of the coalescence of binary systems in circular 
orbit using post-Newtonian (PN) approximations to general relativityta. 

The detection of in-spiral signals is carried out by cross-correlating the detector output 
with a discrete set of filtersB@, usually computed within the so-called restricted post-Newtonian 
approximation Q: PN corrections are taken into account in the phase of the wave-form, whereas 
the amplitude is retained at the lowest Newtonian order. Thus, one discards all multipole 
components except the quadrupole one. Such simplification of the filter structure is believed to 
have negligible effects on the detection performances and is not expected to affect appreciably 
the statistical errors in the estimation of the source parameters. 

Here we investigate the effects of the introduction of PN corrections also to the amplitude 
of the wave-form and discuss some implications for signal detection and parameter estimation. 
For sake of simplicity we will assume negligible spins and, as it is always the case for ground 
based experiments, circular orbits. 

2 The wave-form 

The signal produced at the output of an interferometric detector by a gravitational wave of 
polarization amplitudes h + and h x can be written as 



h(t) = F + h+(t) + F x h x (t) , 



(1) 



where F + and F x are the so-called beam pattern functions of the detectori; they depend on the 
location of the source in the sky (#, <p) and the polarization angle ip. If we consider the in-spiral 
of a binary system of masses mi and ni2, h + and h x readla 
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) 2 / 3 , uj is the system orbital frequency and r the source distance; m = mi + 777-2, 



where x = (mu 
fi = 7771777,2/777,, ?7 = fi/m and M. 



H 3 / 5 m 2 / 5 



mr j3/5 are foe total mass, the reduced mass, the 



symmetric mass ratio and the chirp mass, respectively. 

The lower terms of the PN expansion for the plus and cross polarization are given byt 
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where c = cos 1 and s = sin t; t is the angle between the direction of the source and the orbital 
angular momentum, and & is twice the orbital phase. 

If one considers only the first term in (g), corresponding to the Newtonian one, one gets the 
restricted PN approximation. Our goal is to study the full 2 PN wave-form; here we will present 
some preliminary results, where the 0.5 PN corrections to the amplitude are taken into account, 
keeping however the phase at 2 PN order. Although this is a simplified and, to some extent, 
arbitrary choice of signal, all new features of the wave-form are introduced, in particular more 
information about the masses and the position of the source. 

Considering the amplitude through 0.5-PN order, the GW output at the detector can be 
written as 

h (t) = ^( m ^ F ) 2 / 3 {h^(t) + (mnFfW*)®} , (7) 
where F is the quadrupole gravitational wave frequency, i.e., dQ/dt = 2irF, and 
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The Fourier transform of h(t), calculated using the stationary phase approximation!!, reads: 
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In the restricted PN approximation A = 1. Now it contains two additional contributions 
related to the 0.5 PN corrections to the amplitude. Notice the dependency of these two terms 
on sint and y/Y — 4r/: the departure from the value A = 1 increases as i — > ir/2 and r] — > 0. 

3 Formalism 

We denote by h the "true" GW signal and u(6) the family of templates, as a function of the 
parameter vector 6 = (t c ,(j) c ,\). The signal to noise ratio SNR, for optimal filtering, is defined 
as! 

SNR = s](h\h) , (15) 

where ( | ) denotes the usual inner product. The fraction of SNR obtained by cross-correlating 
a template u{0) with h is given by the ambiguity function 
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which depends on the choice of 0. The maximum of the ambiguity function over the whole 
parameter space is defined as the fitting factor 
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The fitting factor is a measure of how well any chosen family of templates fits the signal h. The 
maximization of the ambiguity function over the extrinsic parameters <f> c and t c , phase and time 
of coalescence, is the so-called match 



M(Ai,A2)= max 
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For the set up of the bank of filters, one sets a minimal match i'0 as the match between signal 
and template in the case where the signal is equidistant from all the nearest templates. 



4 Results 



In the results presented in the following, the signal h is computed according to Eq. (13)-(15), 
whereas the template wave-forms are calculated within the usual restricted PN approximation, 



Table 1: Comparison of the values of the fitting factors FF, and the ratios TZ = W (u\u) / (h\h) for the pair of 
masses 0.1-10 Mq, and 1.4-10 Mq, for different orientations in the sky and polarization angles. 



Orientation-polarization 


t = LN 


n i in ft /r 

U.1-1U Mq 

FF K 


i a in a ,f 
1.4-1U Mq 

FF K 


V0, 0, V 


vr/2 


0.9319 0.9319 


0.9549 0.9560 


= tt/6, = tt/4, ip = vr/4 

= 7T/6, = 7T/4, ?/> = 7T/6 

6» = vr/6, = vr/4, V = 


tt/3 
7r/3 
vr/3 


0.9495 0.9495 
0.9498 0.9498 
0.9516 0.9516 


0.9669 0.9677 
0.9671 0.9681 
0.9683 0.9691 


= tt/6, = tt/4, ^ = vr/4 

# = 0, = 0, -0 = 

9 = vr/6, = vr/4, = vr/6 


tt/4 
tt/4 
tt/4 


0.9662 0.9662 
0.9662 0.9662 
0.9664 0.9664 


0.9780 0.9787 
0.9780 0.9787 
0.9782 0.9789 


6> = vr/6, = tt/4, ip = tt/4 

8 = 0, = 0, V = 

6/ = vr/6, = vr/4, = vr/6 

9 = vr/6, = vr/4, = 


tt/6 
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0.9829 0.9829 
0.9829 0.9829 
0.9829 0.9829 
0.9831 0.9831 


0.9890 0.9895 
0.9890 0.9895 
0.9890 0.9895 
0.9891 0.9896 



corresponding to A = 1 in Eq. (14). The noise curve is the one corresponding to the initial 
LIGO configurationS. 

In table |], we give the fitting factors and the ratios of SNR, TZ = y/ (u\u) / (h\h) , for the pair 
of masses 0.1-10 Mq and 1.4-10 Mq, for different orientation and polarization angles. We note 
that the fitting factor reaches a minimum for i = ir/2 as expected, and monotonically increases 
as l — > or l — > 7t; it depends rather weakly on 0, and 0. 

For a fix position in the sky and i = tt/2, we calculate then the fitting factors for different 
mass pairs. The fitting factor varies from 0.87 to 1.0; FF gets smaller as m increases and/or 7] 
decreases, see table ||. 

It is now interesting to investigate the loss of SNR if one uses a restricted PN bank of 
filters to detect signals that include PN amplitude terms. The discrete mesh of filters is nor- 
mally generated in such a way that any signal in the restricted PN plane produces a match 
M(n(Ai), n(A2)) is always larger than a minimum value, usually set to 0.97. Assuming that the 
"true" GW signal, h, is includes PN corrections to the amplitude, so that it lies outside the 
template space, we want to quantify the match M(/i(Ai), u(A2)) between the PN signal and the 
nearest restricted PN template. 

In table ||, we present the results obtained for a system of masses O.1-1OM and different 



Table 2: Fitting factors and the ratios of SNR for l = 7r/2. 
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Table 3: Fitting factors and match when the GW signal and the template have different parameter values. 

Results for 0.1-10M B . 
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Table 4: Measurement errors at SNR=10. 
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choices of i. For each case we calculate the fitting factor between the signal and the family 
of templates; then, using different parameter values Ai and A2, we calculate the match in the 
restricted PN plane M(u(\\), u(A2)), and, for the same parameters, the match between the 
signal and the template M(/i(Ai), u(A2)). 

What we observe is that the "real match" can be approximated as M(h(Xi), u(\2)) ~ 
FF x M(u(A 1 ),n(A 2 )). 

We turn now attention to the issue of estimating the source parameters. It is important 
to notice that for a waveform computed taking into account PN corrections to the amplitude, 
two additional parameters are involved; our choice corresponds to F x /F + and cos 1. We adopt 
here the standard variance-covariance matrix analysis afij, although it can underestimate the 
statistical errors in the limit of low SNrElIIiI; such well-known problem is beyond the purposes 
of this work and for sake of simplicity we will remain into the usual frame of the computation 
of the Fisher information matrix. 

In table ||, we compare the errors for SNR=10 between the restricted and the "non-restricted" 
PN approximation. The results refer to a system of 1.4-2OM0 and 9 = it/ 6, 4> = vr/4, ip = vr/4, 
and 1 = 7r/3. The use of more accurate wave- forms leads to smaller errors (by roughly a factor 
~ 2) in the determination of the masses; however the information about 1 and F x /F + remain 



Table 5: Correlation coefficients, (a) Restricted 2 PN, (b) 0.5-2PN. 
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very poor. 

In table [| we provide the correlation coefficients for the same parameter choice. We notice 
that the correlation coefficients are also smaller and that the two new parameters have small 
correlations with the other ones. Notice that the amplitude is now correlated with the other 
parameters while it is not the case for the restricted PN case. 
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